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Binary-driven-hypernovae (BdHNe) within the induced gravitational collapse (IGC) paradigm 
have been introduced to explain energetic ( Ei so > 10° 2 erg), long gamma-ray bursts (GRBs) associ¬ 
ated with type Ic supernovae (SNe). The progenitor is a tight binary composed of a carbon-oxygen 
(CO) core and a neutron star (NS) companion, a subclass of the newly proposed “ultra-stripped” 
binaries. The CO-NS short-period orbit causes the NS to accrete appriciable matter from the SN 
ejecta when the CO core collapses, ultimately causing it to collapse to a black hole (BH) and pro¬ 
ducing a GRB. These tight binaries evolve through the SN explosion very differently than compact 
binaries studied in population synthesis calculations. First, the hypercritical accretion onto the 
NS companion alters both the mass and momentum of the binary. Second, because the explosion 
timescale is on par with the orbital period, the mass ejection can not be assumed to be instanta¬ 
neous. This dramatically affects the post-SN fate of the binary. Finally, the bow shock created 
as the accreting NS plows through the SN ejecta transfers angular momentum, braking the orbit. 
These systems remain bound even if a large fraction of the binary mass is lost in the explosion (well 
above the canonical 50% limit), and even large kicks are unlikely to unbind the system. Indeed, 
BdHNe produce a new family of NS-BH binaries unaccounted for in current population synthesis 
analyses and, although they may be rare, the fact that nearly 100% remain bound implies they may 
play an important role in the compact merger rate, important for gravitational waves (GWs) that, 
in turn, can produce a new class of ultrashort GRBs. 
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INTRODUCTION 

Binary massive star systems evolve into a broad set of 
compact binaries from X-ray binaries consisting of stars 
accreting onto either BH or NS companions to the more 
exotic binary compact objects such NS-BH and NS-NS 
binaries. The formation scenarios of these compact bi¬ 
naries typically argue that, after the first SN explosion, 
the compact remnant enters a common envelope phase 
with its companion, tightening the orbit. If the system 
remains bound after the companion star collapses, a NS- 
BH or NS-NS binary is formed. A range of scenarios have 
been invoked HH3j, including exotic scenarios where both 
components expand off the main sequence concurrently, 
causing a single common envelope around two helium 
cores [3J. 

Recently, two independent communities have argued 
for a “new” evolutionary scenario forming these compact 
binaries where, after the collapse of the primary star to a 
NS, the system undergoes a series of mass transfer phases, 
ejecting both the hydrogen and helium shells of the sec¬ 
ondary to produce a binary composed of a massive CO 
core and a NS. When the CO core collapses and produces 
a SN explosion, a compact binary system is formed. In 
the X-ray binary/SN community, these systems are called 
“ultra-stripped” binaries. In the past few years, such sys¬ 


tems have been invoked to both explain the population of 
NS-NS binaries as well as a growing set of low-luminosity 
and/or rapid decay-rate SNe 0[6]. Low-mass ejecta can 
match the observational features of these SNe and ultra- 
stripped binaries without hydrogen and helium layers in 
their pre-SN progenitor produce small cores with such 
low-mass explosions. The rate of these systems are pre¬ 
dicted to be 0.1-1% of the total SN rate [5]. These bina¬ 
ries are extremely tight, and most of the systems studied 
have orbital periods lying between 3000 and 300,000 s. 
Proponents of the ultra-stripped systems argue that this 
scenario dominates the formation of NS-NS binaries and 
that there are virtually no systems that are formed where 
the CO core collapses directly to a BH. 

The IGC scenario for GRBs [ZHS] introduced a sub¬ 
set of extremely short-period CO-NS binaries where the 
ejecta from the exploding CO star accretes onto its NS 
companion, causing the NS, in some cases, to collapse 
to a BH. If ultra-stripped binaries dominate the forma¬ 
tion of NS-NS binaries, this scenario would dominate the 
formation of NS-BH binaries. This collapse to a BH re¬ 
leases energy to drive the GRB emission 0 [TO]. The 
CO core is a requirement to allow the tight orbits needed 
to produce sufficient accretion to cause the NS collapse, 
but it also provides a natural explanation for the fact 
that these GRBs are always associated with type Ic SNe. 
The recently introduced ultra-stripped binaries are a wel- 
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come support for the IGC scenario from the point of view 
of stellar evolution, with the only caveat that IGC pro¬ 
genitors are a small subset of the ultra-stripped binaries 
where the initial orbital separation and CO core mass 
are aligned to produce binaries with orbital periods lying 
in the 100-1000 s range. This requires fine-tuning both 
of the CO star mass and the binary orbit. From an as- 
trophysical point of view the IGC scenario is uniquely 
characterized by the formation of the BH during the ac¬ 
cretion process of the SN ejecta onto the companion NS 
and the associated GRB emission. Since the rate of the 
high-luminosity GRBs (BdHNe) explained through the 
IGC scenario is (1.1 1.3) x 10 -2 Gpc -3 y _1 [IT], and 
0.1-1% of the SN Ibc population could be ultra-stripped 
binaries 0, only 0.005-0.07% of the latter are needed to 
explain the BdHNe population (assuming a SN Ibc rate 
of 2 x 10 4 Gpc -3 y _1 [T2]). 

Studies of ultra-stripped binaries have expanded our 
understanding of stellar radii, confirming these results: 
CO cores with masses below 2 Mq have radii of 1 
4 x 10 9 cm m , in agreement with the assumptions used 
in IGC studies m- Even if some helium remains on the 
stripped core, it will be ejected if it expands to inter¬ 
act with its compact-object companion. These radii are 
sufficiently small to produce the tight orbits required to 
produce the rapid accretion of the ejecta onto the NS 
companion and the formation of the BH. 

In typical systems, most of the binaries become un¬ 
bound during the SN explosion because of the ejected 
mass and momentum imparted (kick) on the newly 
formed compact object in the explosion of the massive 
star. Under the instantaneous explosion assumption, if 
half of the binary system’s mass is lost in the SN explo¬ 
sion, the system is disrupted, forming two single compact 
objects. Although SN kicks may allow some systems to 
remain bound, in general, these kicks unbind even more 
systems. In general, it is believed that the fraction of 
massive binaries that can produce double compact ob¬ 
ject binaries is low: ~0.001 -1% [THS]. 

For ultra-stripped binaries, the fate is very different. 
In these systems, the mass ejected is extremely low and, 
if the SN kick is low, these systems remain bound P3G3- 
In the tighter binaries leading to IGC progenitors, the 
assumption of instantaneous mass ejection is no longer 
valid. We demonstrate in this work that, removing this 
assumption, even with a strong SN kick nearly all of these 
systems will remain bound. In this case, even though 
IGC progenitors are rare, the compact binaries produced 
by these progenitors may dominate the total NS-BH bi¬ 
naries in the Universe, and lead to a new previously un¬ 
accounted family of GRBs. 

We shall describe below the differences between these 
systems and typical massive star binaries, modeling these 
orbits through the SN explosion. We then calculate the 
evolution of these NS-BH binaries via GWs emission up 
to the merger point, and assess their detectability. We 


conclude with a discussion of the additional observational 
predictions of these NS-BH binaries, introducing a new 
class of short GRBs, with specific observational signa¬ 
tures, here referred to as ultrashort GRBs. 

POST-EXPLOSION ORBITS 

The mass ejected during the SN alters the binary orbit, 
causing it to become wider and more eccentric. Assum¬ 
ing that the mass is ejected instantaneously, the post¬ 
explosion semi-major axis is a/ao = (Mq — AM)/(Mq — 
2aoAM/r), where oo and a are the initial and final semi¬ 
major axes respectively, M 0 is the total initial mass of the 
binary system, AM is the change of mass (equal to the 
amount of mass ejected in the SN), and r is the orbital 
separation at the time of explosion II]. For circular or¬ 
bits, like the ones expected from our systems after going 
through a common envelope evolution, we find that the 
system is unbound if it loses half of its mass. But, for 
these close binaries, a number of additional effects can 
alter the fate of the binary. 

The time it takes for the ejecta to flow past a compan¬ 
ion in a SN is roughly 10-1000 s. These explosions follow 
a so-called homologous velocity profile where the velocity 
is proportional to the position. Although the shock front 
is moving above 10,000 km s _1 , the denser,lower-velocity 
ejecta can be moving at below lOOOkms -1 . Our esti¬ 
mates are based on simulated supernova explosions IIP] , 
The broad range of times arises because the SN ejecta 
velocities varies from 100-10,000 km s _1 . The accretion 
peaks as the slow-moving (inner) ejecta flows past the 
NS companion. Note that the initial SN explosion in 
this case is not a hypernova. The observed “hypernova” 
is actually produced when the GRB from the BH col¬ 
lapse sweeps up this SN (and circumstellar) material [15] . 
For normal binaries, this time is a very small fraction of 
the orbital period and the “instantaneous” assumption 
is perfectly valid. However, in the close binary systems 
considered here, the orbital period ranges from only 100- 
1000 s, and the mass loss from the SN explosion can no 
longer be assumed to be instantaneous. 

This has already been pointed out in |}16j where it was 
shown that in BdHNe the accretion process is fast and 
massive enough to produce the BH formation in a time- 
interval as short as the orbital period. We here deepen 
this analysis to study the effect of the SN explosion in 
such a scenario with a specific example, for which we 
have produced an orbit code using a simple staggered 
leapfrog integration (see E3 for details of this integra¬ 
tion method). We have tested both stability (by mod¬ 
eling many orbits) and convergence (decreasing the time 
step by 2 orders of magnitude confirming identical re¬ 
sults) . We also reproduce the results of the instantaneous 
limit. From figure [I] as the ejecta timescale becomes 
just a fraction of the orbital timescale, the fate of the 
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FIG. 1. Semi-major axis versus explosion time for 3 differ¬ 
ent mass ejecta scenarios: 3.5 Mg (solid), 5.0 Mg (dotted), 
8.0 Mg (dashed). The CO core collapse to form a 1.5 Mg 
NS (its initial mass is the ejecta mass plus the NS mass), 
and the companion NS has a mass of 2.0 Mg. If the explo¬ 
sion were instantaneous, all of our systems with ejecta masses 
above 3.5 Mg would be unbound. For explosion times above 
1.2 times the orbital time, not only are the systems bound, 
but the final orbital semi-major axis is less than 10 times the 
initial separation. 

post-explosion binary can be radically altered. For these 
models, we assumed very close binaries with an initial 
orbital separation of 7 x 10 9 cm in circular orbits (such 
close binaries are only formed through a common enve¬ 
lope phase which circularizes the orbit). With CO core 
radii of 1-4 x 10 9 cm m, such a separation is small, 
but achievable. We assume the binary consists of a CO 
core and a 2.0 M 0 NS companion. When the CO core 
collapses, it forms a 1.5 M 0 NS, ejecting the rest of the 
core. We then vary the ejecta mass and time required 
for most of the ejected matter to move out of the binary. 
Note that even if 70% of the mass is lost from the sys¬ 
tem (the 8 Mq ejecta case), the system remains bound 
as long as the explosion time is just above the orbital 
time (T or bit = 180 s) with semi-major axes of less than 
10 11 cm. 

The short orbits (on ejecta timescales) are not the only 
feature of these binaries that alters the post-explosion 
orbit. The NS companion accretes both matter and mo¬ 
mentum from the SN ejecta, reducing the mass lost from 
the system with respect to typical binaries with larger 
orbital separations and much less accretion. In addition, 
as with common envelope scenarios, the bow shock pro¬ 
duced by the accreting NS transfers orbital energy into 
the SN ejecta. In figure [2] we show the final orbital sep¬ 
aration of our same three binaries, including the effects 
of mass accretion (we assume 0.5 Mq is accreted with 
the momentum of the SN material) and orbit coupling 
(30% of the orbital velocity is lost per orbit). With these 
effects, not only do the systems remain bound even for 
explosion times greater than 1/2 the orbital period but, 


FIG. 2. Semi-major axis versus explosion time for the same 
3 binary systems as in figure [l] including mass accretion and 
momentum effects. Including these effects, all systems with 
explosion times above 0.7 times the orbital time are bound 
and the final separations are on par with the initial separa¬ 
tions. 

if the explosion time is long, the final semi-major axis 
can be on par with the initial orbital separation. 

The tight separation of these binaries facilitates tidal 
locking and the angular momentum axis of the CO core 
will be aligned with the orbital angular momentum. For 
many of the kick mechanisms in the literature, the kick is 
often aligned with the rotation axis. For example, both 
in neutrino-driven mechanisms USES and asymmetric 
explosions driven by convection pi2 [2(1, 21) the kick is 
aligned with the rotation axis. However, it is still possi¬ 
ble to have some misalignment leading to some eccentric¬ 
ity and “tumbling” of the system with specific signatures 
in the light curve following the prompt emission of the 
GRB. Hence, we here consider both kicks aligned with 
the rotation (and hence orbital) axis as well as random 
kicks. If the kick is aligned with the orbital plane, the sys¬ 
tem can remain bound even with kick velocities as high 
as 1000 km s -1 . However, if the kick is in the same direc¬ 
tion as the star is moving, the systems can be disrupted 
if the kick is above 500-700 km s -1 if the accretion phase 
is longer than an orbital period. 

The tight compact binaries produced in these explo¬ 
sions will emit GW emission, ultimately causing the sys¬ 
tem to merge. For typical massive star binaries, the 
merger time is many Myr. For BdHNe, the merger time 
is typically 10,000 y, or less (figure [3]). 


GRAVITATIONAL WAVES FROM THE NS-BH 
BINARY 

To better understand the GW signal from these merg¬ 
ers, we study the evolution of the orbital binding en¬ 
ergy Eb up to the merger following the effective one-body 
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FIG. 3. Merger time due to GW emission as a function of 
explosion time for the same 3 binary as in figure [l] includ¬ 
ing mass accretion and momentum effects. Beyond a critical 
explosion time (0.1-0.6 Torbit depending on the system), the 
merger time is less than roughly 10,000 y. For most of our 
systems, the explosion time is above this limit and we expect 
most of these systems to merge quickly. 


(EOB) formalism I221 - 1T5| up to the 4th Post-Newtonian 
approximation (see Refs. EH 121 ] and references therein). 
We adopt here Mns = 1-5 M© and Mbh = 2.67 M 0 
[25] . the latter corresponding to the critical mass M cr i t 
of a non-rotating NS obeying the nuclear NL3 equation 
of state (EOS). Uncertainties in the EOS at supranuclear 
densities lead to a variety of NS mass-radius relations and 
consequently to different values of M cr j t , hence of Mbh- 
Both rotation [29] and different binary parameters may 
lead to different amounts of angular momentum trans¬ 
ferred to the NS, affecting its mass [30 . 

In order to assess the detectability of the GW emission 
by advanced LIGO (aLIGO), we compute the signal-to- 
noise ratio (SNR), averaged over all sky locations and 
binary orientations, (SNR), generated by the NS-BH 
spiraling-in binary up to the merger point m Follow¬ 
ing [SB, we adopt as a threshold for the aLIGO detec¬ 
tion (SNR) = 8 in a single detector, which implies a GW 
horizon distance for these NS-BH binaries, which have 
a chirp mass M ch = (M B hM N s) 3 / 5 /(M B h + M NS ) 1/5 ~ 
1.73 M©, d,L « 335.4 Mpc, or zk 0.075, using the maxi¬ 
mum possible sensitive reachable by 2022. No BdHN has 
been up to now detected with such a low redshift. Fig¬ 
ure [4] shows, for two sources shown to be consistent with 
the BdHN picture (GRB 130427A with z = 0.34 [T5j and 
GRB 061121 with z = 1.31 [22]), the GW source ampli¬ 
tude spectral density, = 2\h(f d )\y/fd = M/d)/\/7d, 
together with the one-sided ASD of the aLIGO noise, 
VSn(fd)- Here h c (fd) and h(fd) are the characteristic 
strain and the Fourier transform of the signal and fd the 
frequency of the GWs at the detector. For these sources, 
(SNR) ss 1.75 and 0.45, respectively. For an optimally 
located and polarized source, the SNR could increase by 


up to a factor ss 2.26, which implies that SNR=8 could be 
obtained for a source as far as d,L ~ 2.26 x 335.4 Mpc ss 
758 Mpc, or z « 0.160. Furthermore, the SNR scales as 
A4^ 6 , so it increases e.g. with larger BH masses. For 
rotating NS with the NL3 EOS, the maximum value of 
M cr it is ~ 3.4 M© :22], which would increase the SNR 
only by « 1.1. For this largest BH mass, the GW horizon 
becomes ~ 1.1 x 758 Mpc ~ 834 Mpc, or z w 0.174. 
This largest possible GW horizon implies an upper limit 
of ~ 0.03 detections per year, adopting a BdHN rate of 
1.2 x 10 -2 Gpc -3 y _1 [IT] . 



FIG. 4. ASD of the spiraling-in phase up to the merger, 
= 2|ft(/d)|v7a = h c {f d )/y/Td, of the NS-BH binaries 
produced by two BdHNe, GRB 130427A at redshift z = 0.34 
and GRB 061121 at z = 1.31, compared with the noise ASD of 
aLIGO, vSW(/ti). We indicate the estimated SNR for these 
two sources and show the case of the NS-BH binary which 
would generate a positive detection with SNR=8. The binary 
dynamics is simulated via the EOB formalism up to the 4th 
Post-Newtonian approximation. 


CONCLUSION 

The evolutionary scenario for BdHNe requires much 
tighter binaries than typically studied in the literature of 
ultra-stripped binaries and this produces unique features 
in the end-fate of these systems. The progenitor of this 
GRB engine begins with two massive stars, in contrast to 
the one based on a massive core collapsing to a BH [33]. 
A tight binary is produced after a succession of common 
envelope phases, producing a CO core near Roche-Lobe 
overflow orbiting a NS, a subset of the ultra-stripped bi¬ 
naries j5, 6] [34] 35]. Since 0.1-1% of the total SN Ibc 
are expected to be ultra-stripped binaries [5], we esti¬ 
mate that only 0.005-0.07% of the latter are needed to 
explain the observed population of BdHNe. The fate of 
such systems evolves very differently than the standard 
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picture. The NS can accrete appreciable material in the 
SN explosion and this accretion causes it to collapse to 
a BH and form a GRB. However, the tight binary inval¬ 
idates many of the assumptions about orbital evolution 
in the SN. The SN explosion does not pass “instanta¬ 
neously” across the NS, and correcting this assumption 
alone drastically alters the binary fate. Including the in¬ 
teraction of the orbit and the ejecta further exacerbates 
these differences, causing these NS-BH to be very differ¬ 
ent than the systems prediction in standard population 
synthesis models. 

First and foremost, the fraction of the BdHNe that re¬ 
main bound after the SN explosion is nearly 100% even 
with large ~ 500-1000 km s -1 kicks imparted during 
the SN explosion instead of the <1% in standard scenar¬ 
ios. This means that even if BdHNe are rare, they may 
dominate the fraction of NS-BH binaries in the Universe. 
In addition, the merger timescales for these systems are 
typically <10,000 y, producing a set of rapidly-merging 
binaries. In view of such a short lifetime due to GW 
emission, the current number of such events is likely to 
be comparable with the original rate of long GRBs pro¬ 
duced by BdHNe following the IGC paradigm. Because 
of this rapid merger, the systems are unlikely to travel 
that far from the site of the SN explosion that formed 
the GRB. Even with large kicks, we expect these bina¬ 
ries to merge within 10 pc of the BdHNe and we expect 
the merger to occur within the radius swept clean by 
the BdHN, giving a characteristic imprint in the GRB 
emission. In view of the expected paucity of the bary- 
onic contamination around the merger site, it is expected 
that the characteristic prompt radiation emission time of 
the GRB produced by these sources be dominated by the 
general relativistic timescale of the BH, GM/c 3 ~ 10 -4 - 
10~ 5 s, which justifies the attribution of the name of 
ultrashort GRBs to this new family of events. 

Another observational feature of these binaries is that 
the BHs from these systems are low mass: ~3-4M 0 , of 
the order of the critical mass of rotating NSs 1221 EDI, 
instead of the 5-10 M 0 produced by standard scenarios. 
However, further accretion of mass and angular momen¬ 
tum from material kept bound into the system after 
the BdHN process might lead the BH to larger masses 
and to approach maximal rotation [SU]. Although the 
NS in this NS-BH binary should be rapidly rotating, 
producing pulsed emission, the short timescale between 
formation and merger means that it will be difficult to 
observe such systems through steady pulsed emission. 
However, if these systems make up a sizable fraction 
of the NS-BH population, they could be detecting by 
their GW signal. Although it is difficult to get the 
exact component masses from aLIGO, evidence [36] , or 
the lack thereof, for binaries with low-mass BHs could 
support, or limit the rate of, this scenario. 
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